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Abstract 

We consider a particle system studied by E. Brunet and B. Derrida [9], which 
evolves according to a branching mechanism with selection of the fittest keeping the 
population size fixed and equal to N. The particles remain grouped and move like a 
travelling front driven by a random noise with a deterministic speed. Because of its 
mean-field structure, the model can be further analysed as N —» oo. We focus on 
the case where the noise lies in the max-domain of attraction of the Weibull extreme 
value distribution and show that under mild conditions the correction to the speed 
has universal features depending on the tail probabilities. 
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1 Introduction and main result 

We consider a model of front propagation introduced by Brunet and Derrida [9j. A constant 
number N of particles evolve on the real line in discrete time. Let A'i(0),..., X/v(0) be the 
particles initial positions. With (s); 1 < i,j < N, s > 1} an i.i.d. family of r.r.v.s, the 
positions evolve as 

Xfit + 1) := max {Xj(t) + £jfit + 1)}. (1.1) 

For € L 1 , it is proved in [10] that the limits 

1 1 

lim — max Xfit) = lim — min Xfit) = vm(£) 

i->oo t l<i<N W i-xx> t l<i<N W 

exist a.s. with u/v(£) a real constant depending on the law of £. The limit u/v(£) is called the 
speed of the IV-particle system, and we study here its asymptotic for large N. 
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The model can be viewed as the combination of a branching step with a fixed number N of 
offspring per individual each one being subject to mutation given by the addition of a fresh r.v. 
£, and a selection step where only the N “rightmost” among these N 2 offspring are kept. One 
major question of the held is to understand what really determines the motion and derive the 
universality properties of such models. In particular, the speed of propagation depends on the 
parameter tuning how stringent the selection is, and one is interested in the corrections with 
respect to the speed of the model without selection. The definition of “rightmost” used in the 
present paper, see (EH), is somewhat specific, and is different from the traditional choice of M- 
branching random walks urn when all newborn individuals are simultaneously compared. For 
the latter choice we mention [3j and ini, and also m for the continuous case. A dual problem 
is the survival of the branching population killed by a moving obstacle, e.g. a line {!]. In its 
general form, the model relates to propagations of pulled fronts, when the motion is determined 
by the leading edge m- Archetypes of pulled fronts are branching random walks or branching 
Brownian motions, described by the Kolmogorov-Petrovskii-Piscunov equation. There, and in 
contrast to the present case, one looks for the second order correction in time of the rightmost 
position to its leading order [6]. Though effective equations in the continuum are available to 
describe front dynamics, the process here is intrinsically random and discrete, adding interest 
to its understanding. We note from [16J that asynchronous dynamics leads to free boundary 
problems. 

Already mentioned in Em, the model ( 11.111 was taken up in [9j and studied in the case of 
Gumbel distribution for £, which leads to an exact solution for fixed N, and results have been 
extended in a perturbative picture provided that £ has an exponential upper tail m- In the 
present paper, we consider perturbations of the Weibull distribution, including bounded £’s with 
a polynomial density close to its maximal value. 

From a different perspective, our model can be interpreted as a first passage percolation. By 
a simple induction argument, one obtains from m the formula 

t 

Xi(t ) = max |A J0 (0) + ^ £j s _ Ua 0); 1 < j. s < N, Vs = 0,... ,t-l and j t = i}, (1.2) 

S=1 

which yields a path representation of the interacting particle system. We interpret now —£ji(t+ 1) 
as the passage time on the oriented edge from (J,t) to (i,t + 1): As (11.21) shows, the negative 
of Xi(t) is the passage time from the line t = 0 to the point in a model of first passage 

percolation on the vertex set {1,... , N} x N, and ujv is the so-called time constant of the model. 
Here the graph is oriented (t-coordinate increases by one unit at each step of the path), though 
on the transverse direction jumps are allowed between all pairs of sites j. i (1 < i, j< N). Since 
the graph is complete in the transverse direction, the model is of mean-field type. For general 
percolation models the value of the time constant is not available, but in the present case the 
mean-field feature allows us to determine the time-constant up to first order in the limit of large 
graphs. In the particular case of exponential passage times, the first formula in ( 11 . 81 ) below is 
in force. 

To describe our framework, denote by A (u) the logarithmic generating function of £jj, 

A(u) := lnE[exp(u£ij)], 

and let T>\ := {u £ M; A (u) < oo} be its domain. In this paper, we will assume that the following 
hypothesis hold: 
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Figure 1: Cramer transform I g and Vn- 


(HI) 0 G T>^ (the interior of T>\). In particular, has finite moments of all orders. 

(H2) For every N G N there exists a un G fl [0, oo) such that 

unA'(un) — A(iqv) = In N. (1-3) 


The function uA'(u) — A(u) is increasing on n [0,oo), hence un is unique. Under these 
hypothesis the number 

v N ■■= A'(u n ) (1.4) 

is well defined. If I^(v) is the Cramer transform of £ij 

I^{y) := sup { vx — A(x)}, 


then vn is determined by I^(vn) = lnN,VN > E[^], see Figure HI and it holds I'c(vn) = un. 

In Section [3l we show that vn is an upper bound for the velocity vn(0 of the IV-particle 
system. To obtain a lower bound to w/v(£), we do some additional assumptions on ^ and 
focus on a more restrictive class of distributions. Denoting by F(-) the common probability 
distribution function 

F(x) := P(£y < x ), 

we will further assume that F(-) belongs to the max-domain of attraction of the Type III extreme 
value distribution with probability distribution function 'F Q (-) given by 


( exp (—|x|“) if x < 0 
\ 1 if x > 0, 


(1.5) 


for some a > 0. This law is sometimes called reverse-Weibull (see e.g. Chapter 1 in M) , or 
Weibull for short, and it is the law of —£ l / a with £ an exponential variable with mean 1. It is 
well known that F(-) belongs to the domain of attraction of 4' cl (-) if and only F(-) has a finite 
right-end 

x^ := sup{x G M; F(x) < 1} < oo, 
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and for each x > 0 


lim 


1 — F[x£ — hx) 


= x 


h^>o+ 1 — F (xg — h) 

see, for example, Proposition 1.13 in Section 1.3 of [22]. In this case, let 

ajsr := x£ — inf{.x: F(x) > 1 — 1 /IV}, (1.6) 

then F n (x£ + xcln) —>• ^(x) as N —»• oo and 

lim Nil — Fixe — ojv)) = 1 

N-> oo v 5 ' 

The main result of this paper is the following theorem concerning the speed of the N particle 
system. 


Theorem 1.1. Assume that (HI), (H2) hold, and that ^ij belongs to the domain of attraction 
of the extreme value distribution A> a , for some a > 0. Let 


a 


Cry • 


a)a) 


(1.7) 


where T(-) is Euler’s gamma function and e = 2.718 ... is the Napier’s constant. Then, the speed 
vn{() of the N-particle system satisfies 

v n(£) = X£ — c a aN + o (ajv) as N —> oo, 

where a at is given by 11.6]) . 


Theorem 11.11 gives a partial answer to the problem raised by E. Brunet and B. Derrida 
[9], and we obtain the finite-size corrections to the speed for a large class of distributions that 
are bounded from above. Our result comprises, for example, the negative of the exponential 
distribution and the uniform distribution, for which the correction to the speed are 

v n {~£) ~ —^7 and ( v N (U) - 1) ~ —^7 as N ->• 00 , (1.8) 

eiv eN 

respectively. In the above formulas, means that the ratio of the sides approaches to one as 
IV —>• 00. 

Warm-up calculations: Let us explain how to determine the order of magnitude of the 
correction from elementary considerations. Assume in this paragraph that x^ = 0. On the one 
hand, we can bound from below our IV-particle system with a single particle following the leader, 
i.e., the random walk with jumps law given by maxj<Ar £ii, resulting with a lower bound for u/v(£) 
of order ajy. On the other hand, a naive upper bound is given by the random walk with jumps 
maxjjcjv fji, which leads to a different order 0{a^ 2 ) of the correction for the maximum is over 
N 2 variables this time. One can improve the upper bound by using the first moment method 
of Section [3] leading to the same order O(oat) as the lower bound. However, the multiplicative 
factors do not match, and some deeper understanding and improvement of the lower bound is 
needed. This is what we implement in Section [I] using a comparison with a branching random 
walk with selection. 

Organization of the paper: In Section [2] we present some point processes and branching 
random walks related to our model and we sumarize the necessary results for our purpose. We 
prove the upper bound for the speed in Section [3] by a first moment estimate, and the lower 
bound in Section [4] by coupling, the two bounds resulting in Theorem 11.11 
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2 Point processes and branching random walks 


In this section, we introduce different processes entering the analysis of the particle system m- 

2.1 Point measures on R 

It is convenient to represent populations of particles by point measures on R. Given a vector 
x € R n with coordinates x\,..., x n , one can associate the point measure 

n 

* := £<w 

We use the notation to represent the set of all simple point measures on R, which are locally 
finite and have a maximum.Throughout this paper, a point process is any random variable C 
taking values on Mb- 

Conversely, an element v € M.b can be described as a sequence v = (z/j) i=1 2) ... (possibly 
finite) such that 

v\>v 2 >... 

We denote by max(u) = v\ the maximum of the support of v, and by \v\ = u(R) < oo the number 
of points in v. If two point measures v and fi have the same number of points \v\ = \fi\ = I\ , 
we can define the distance 

\\v~n\\= sup {\vi - im\}. (2.1) 

l<i<K 

We use the notation to denote the usual stochastic ordering 

v -< /j, if and only if v[x, +oo) < n\x, +oo); Vx € R, 

and we will say that V bounds /i from below”, in this case \v\ < \n\. If we represent v and fi as 
an ordered sequence of points, then v -< n implies that 

Vi < ^ for every i < \v\. 

With a slight abuse of notation we will say that the vector x € R n bounds y € R m from below 
and denote “x -< y" if the point measures ac, y associated to x and y respectively satisfy ?c ■< y- 

2.2 Poisson point processes on ] - oo,0] 

In this section, we present some elementary facts concerning Poisson Point Process 


& = { > ^2 > • • • } C R_ , 

with intensity measure \z\^Cdz on R_; we use the abbreviation PPP and assume that C > 
0 ,f3 > —1. For K > 1, the point process 

^ {K) ■■= (&i)i< K ( 2 - 2 ) 

consisting in the K largest points of & will play an important role in the next sections. 
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For C a random point measure on M_, we denote by ip{u \ C) its logarithmic moments 
generating function 

V >(u | £) := InE 

We can easily compute the logarithmic generating function of the PPP and of its A-truncation. 



Lemma 2.1. For f3 > — 1 ,C > 0, let F? be the Poisson point process on (—oo,0] with intensity 
measure qi(dz) = \z\^Cdz, and its largest K points. For u > 0 we have 


and 



— OO 


T(1 + /3)C 

yl+P ’ 



J — OO 


-T(2 + /3)C 
u 2+p 


r k 


lim E 

oo 


yy e u&>i 

_ i=\ 


T(1 + /3)C 

u 1+ P ’ 


lim E 

if—^oo 


K 

£■ 

,i =1 




-r (2 + /3)C 


Proof. The first claim is obtained by the Campbell formula (see Chapter 9 in 11 ~>! ) and the 
second claim is obtained by monotone convergence. □ 


Corollary 2.2. Under the assumptions of Lemma \2.1l the sequence ip(u \ converges 

uniformly on the compacts o/M + to ip (u | as K —>• oo. In particular, if uk > 0 is such that 


iP{uk | ^ K ' > ) = iP'(uk | &W)u K 


then 

(“* 1 ~ ( crdVffl ) 7 “ ' (2 ' 3) 

Proof. The compact convergence is a direct consequence of the pointwise convergence together 
with the monotonicity of ip(u \ and the continuity of ^(it | in u (Dini’s theorem). 

Let 

u 00 = e(r(l + /3)C) 1/1+/3 

then from the first part of Lemma 12.11 we have that 

ipiuoo | &) = ip'iuoo | &) Uoo and ip'(uoo \ &) = - 1 ^ (c , r(l + /3)) 

By point 2 of Lemma [2.11 uk —>• u 00 and the second claim follows from the uniform convergence. 

□ 


2.3 Branching Random Walks 

Branching Random Walks (BRW for short) have been extensively studied in the past years, see 
the seminal paper [T] for a general literature and important results on the subject. In this paper, 
we focus on BRW defined as follows. Let £ be a point process on M, which defines the positions 
of particles and the reproductive law of the underlying Galton-Watson tree. The process starts 
with one particle located at 0. At each time step t —>• t + 1, the particles of generation t die 
and give birth to independent copies of C, translated by their position. We use the notation 
BRW(£) to denote a BRW defined as above. 
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Let T be the Galton-Watson tree obtained by the genealogical tree of the process, thus, its 
offspring distribution is \C\. To each point (or individual) of the BRW(£) one can associate 
a unique vertex w E T. Let e € T be the root of the Galton-Watson tree, then for a vertex 
w E T, let \e,w\ denote the shortest path connecting e with w, and |u;| the length of this path. 
We will sometimes write its points |e,ic] = (e, w\,... ,Wk) with i = |rcj| and = w. It is a 
standard property of Galton-Watson trees that the process starting from a vertex w E T is also 
a Galton-Watson tree with the same distribution. For two vertices w and w' in T we denote by 
ww' the vertex in T in generation |u>| + \w'\ obtained by concatenation. 

We also denote by r](w) the positions of the individual w£l, and by y(t) the point measure 
associated to the BRW(£) 

y(f) 'y ] < ^{^(«')} • 

w^T\\w\=t 

Finally, an infinite ray [e, Woo] := {e,w±,W 2 , ■ ■ ■} C T is an infinite collection of vertices (or 
infinite path), such that Wi is the parent of It represents a family branch in the BRW that 

has not extinguished, and is parametrized by an element E dT of the topological boundary 
<9T of the tree. 

Under mild conditions on £, the asymptotic behaviour of max ( y(t )) is known [HE], that 
we recall now. Assume that for some a > 0, 


E[|£| 1+a ] < oo, (2.4) 

a condition which can be weakened [l], but in this paper we will always have \C\ = K a constant, 
which trivially implies (12.41) . We also assume that the logarithmic generating function for the 
branching random walk 

ip(u | £) := InE 

is finite in a neighbourhood of u = 0 and that there exists a u* = u(C) > 0 for which 



i/>(u* | c) =«V(«* I £). 


( 2 . 6 ) 


If (12.41 - 12.61) hold, there exists a constant 7 (£) depending only on £ such that 


lim - max (y(t)) = y(£) 

t ~*OO t 


a.s. 


(2.7) 


Moreover, the asymptotic speed 7 (E) is explicit and given by 7 (£) = il>' ( u* \ £), see [H 21 Fl7| 
for a rigorous proof and more results on the subject. 

The theorem below, proved by Gantert, Hu and Shi im gives the precise decay for the 
probability that there exists an infinite ray in the BRW that always stays close to y(£). 

Theorem 2.3 ([17]). Let C be a point process satisfying \2.f \ - 1 2. hi) and (r](w);w E T) be the 
BRW(C). Given 5 > 0, denote by p(oo,5) the probability that there exists an infinite ray in the 
branching random walk that always lies above the line of slope 7 (E) — 5. 


p{ 00 ,5) := E dT : rj(w t ) > ( 7 (£) - 5)t, Vw t E [e, , 


where wt E [e, Woo] is the vertex in generation t. Then, as 5 \ 0 

( ( \x(£) + °(i) n 1/2N 

p{ 00 , 0 ) ~ exp — 


where x(£) = u*if" (u* \ £) for u* given by \2.6\) . 
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2.4 Models with selection and the M-BRW 


Recently, some models of evolving particle systems under the effect of selection have been studied 
[3 ns ns mum ns- The selection creates correlation between individuals of same generation 
and additional dependence in the whole process. 

Berard and Gouere [3] focused on the binary Branching Random Walk with selection of the 
M rightmost individuals (the M-BRW). It consists in a BRW, subject to the effects of selection, 
defined by the point process 

L- Sp 1 -j- d p 2 , 

where pi are i.i.d. As soon as the population size exceeds M, we only keep the M rightmost 
individuals and eliminate the others. If at time zero the number of particles is already M, the 
population size is kept constant. 

Denote by z/m(I) the point process generated by this M-BRW. Berard and Gouere [3] show 
that the support of yM{t) has a diameter of order In M. They also prove that under some 
assumptions on the exponential moment of pi there exists a constant 7 m(£) such that 

lim t~ l min (i/m(£)) = hm t -1 max (j/m(^)) = 7 m(£) a.s. 

t —»oo v ' t—Kx> v 

The existence of 7 m(£) is obtained by the Kingman’s sub-additive ergodic theorem, and by 
monotonicity arguments one can prove that 7 m(£) converges as M —> 00 . The striking result is 
that Berard and Gouere computed the asymptotic limit and the rate of convergence: 

7 m(£) = 7 (£) + x(£)(lnM ) -2 + o((lnM)~ 2 ) as M -> 00 , 

where 7 (£) is the asymptotic speed for the BRW(£), and ,\(£) is from Theorem 12.31 
One can easily define more general M-BRW, let £ be a point process, for which 

\C\ > 1 a.s. (2.8) 

so the process does not extinguish. Under similar assumptions on the exponential moments of 
£, one can prove that the cloud of M points does not spread and that it travels at a certain 
speed 7 m (£), see [19]. 

In Section [I] we show that under the hypothesis of the Theorem ll.il the IV-particle system 
(11.11) can be bounded from below by a family of M-BRW indexed by N. We will then adapt the 
arguments in [3] to obtain a uniform lower bound for the speeds of the BRWs. 

2.5 Elementary properties of Brunet-Derrida’s IV-particle system 

In this section, we present some elementary properties of the IV-particle system. We also intro¬ 
duce some notations that will be useful in the forthcoming sections. Most of these properties 
have been rigorously proved in [ID], therefore we will simply outline the main ideas. 

It will be convenient to consider the process X*(t) obtained by ordering the components of 
X(t) at each time t. Denote by 

X (1) (t) > X (2) (f) > ... > X (N \t) 

the components of X*(t). Let, also a = a(t) be the random permutation of {1,...,IV} such 
that 

X^(t) = X ai{t) (t). 
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Such a ranking permutation is unique up to ties, which we break in the order of the original 
labels. Consider the random variable 

T := inf jf > 1; ^(i- i),i(t) = Vi = 1... ivj. 

Then, T is a stopping time for the filtration 

T t := a{Xj(0), ^j(s); 1 < i,j < N and s < t}. (2.9) 

It has a geometric distribution with parameter not smaller than (1/N) n . Moreover, in generation 
T the position of each particle X t (T) is determined by the position of the leader X ai in generation 
T — 1. We define the process seen from the leading edge 

:=Xi(t)-X ai{t) (t). 

It is Markov process on M. N , which is irreducible, aperiodic and Harris recurrent (due to the 
renewal structure), thus there exists a unique stationary measure ^ r, and for any starting point 
X (0) = x the law of X(t ) converges in total variation distance, 

distT.v. [Jf(X(t) | X(0) = x),7 r) —>• 0, as t —>• oo. 

It proves, in particular, that the cloud of N points remains grouped as t —>• oo. Moreover, by 
the renewal and ergodic theorems, the limit 

1 1 

vn(Q = lim - max {Xi(t)} = lim - nnnjX^t)} 

t —^oo z t —^oo z l<i<N 

exists a.s., see Section 2 in m for a rigorous proof and more details. 


3 Upper bound for the velocity 


In this section, we show that if Cj satisfies the hypothesis |(H1)[ |(H2)[ then 

vn(Q < vn = A'(un), 


where un > 0 is the unique positive solution of uA'(u) — A (u) = In IV. The idea is to use the 
so-called first moment method to bound the probability 


( max |JVj(t)) > tA 1 (un) 
\l<i<N 1 J 


A first and simple observation is that the initial position of the particles does not change 
the speed of the Wparticle system. Hence, we may assume without loss of generality that all 
N particles start at zero. Using the representation (11.21) one gets 


max 

1 <i<N 


W*)} 


max I ^ (s); 1 <j s < X , Vs 

l 8=1 
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By the union bound and Chernoff bound we obtain, for v > vn = A '(un) (which is larger than E[£] 
for A sufficiently large), 


(i<?;< X y - tv ) = P (^do,ji, • • • ’it : - ty ) 


<N t+1 P 

\s=l / 

< N t+1 exp ( — , 


(3.1) 


for all A € N. S ince | (H f)| and [(H2)| hold, I%(v) exists and I$(v) > InA. As a consequence, (13.11) 
has a geometrical decay as t —> oo, which implies, by Borel-Cantelli lemma, 


^ rna x^Xi(t) > tv for infinitely many t € N ) =0, 


hence, lirn sup^^ t 1 riiax{X ? ;(t)} < v a.s. for every v > vn, finally yielding 

limsup- max Xjft) < vn a.s. (3-2) 

t—»oo t l<i<N 


We formalize this result in a proposition. 


Proposition 3.1. Assume that (HI), (H2) hold. Let un > 0 such that unA'(un) ~ A(itjv) = 


In A" and vn = A '(un). Then, for every N € N, 


vn(0 < v N . 


Hence the next step is to study the asymptotic of that we start with the case x% = 0. 

Proposition 3.2. Assume that the hypothesis of Theorem ] 1.1\ hold with x £ = 0. Let itjy > 0 be 
the unique solution of unA'(un) — A (un) = In A, c a be given by \1. 7| ) and ajy by M.h]) . Then, 
as A —» oo 

A'( un ) = ~c a aN + o (ajv), 
which implies that limsup^y -x x > a N 1 VN(£) < — c a . 

Proof. By definition of A(-), we have that 


E [e UiV ] 



In A. 


Note that un —t oo as N —>• oo, indeed it is a direct consequence of the monotonicity and 
continuity of uA'(u) — A (u). Hence, the asymptotic behaviour of the Laplace transform of 
in un is determined by its behaviour in a neighbourhood of zero. Since (ij is in the domain 
of attraction of the function 1 — F{—x) : R + —» M + is a-regularly varying at zero. By 

Karamata’s representation (see Chapter 0 in m) 


1 - F(-x) = Pfej > -x) 


x a c(x x )exp 




x > 0, 
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where c(-) and e(-) are positive functions such that c(t) —>• c > 0 and e(t) —>• 0 as i —>• oo. As a 
consequence, given e > 0 , one can find a u £ > 0 such that for 0 < u < u e 

1 — F(—u) > (c — e)n a . 


Now, we compute the Laplace transform of in un 


E[e UN ^]=¥^ ijUN >-l) I e 

= (1 ~F (-u* 1 )) 


_ z lP(^jUAr > -z) 
0 E(^ijUN > —l) 
V“iv 


dz 


dz + 


■■■ dz 


'y/UN 


We analyse each integral separately. For N sufficiently large < u £ , hence 


-z^{^3 Un — ~ z ) • - u ' 


dz < 


'yfim ^{£,ijUN > -l) ( c_£ ) 


e Z dz, 


/UN 


which converges to zero as N —> oo. Take L > 0, and assume that y/uN > L , then 


-vAW _^P(^-tx JV > -z) 


dz = 


/ 0 > — l) JO 

Using dominated and monotone convergence we obtain that 


dz + 


ry/UN 


■ ■ ■ dz. 


lim lim 


^ e _M ijUN >-z) d \ =T{a + 1) 


irrii xxxxx i f V/ / \ ' 

L—>oo N—hx> yj o P(£yltAT > — l) 

Finally, we prove that the integral from L to ^/njv vanishes as N —>• oo and L —>• oo (in this 
order). For L > 1 and L < z < yTt/v, Karamata’s representation yields 


a c(z 1 1 /jv) 

= z - 7 - 7 — exp 


r u N 


e(t) 


E(£,ijUN > — 1 ) 


c(u N ) 


dt 


IUn /z 


Taking N sufficiently large so e(t) < £ and |c — c(t)\ < £ for every t > Jun 


E(^ijUN > z) < (c + g) ~ £ 
E(CijU N > -1) _ ~ (c - e) 


which yields the upper bound 


> ~z)_ dz < (c + e) 


/; 


_ 2 z a+£ dz. 


A P(£y«JV > —l) ( c ~ e ) Jl 

The right-hand side of this inequality decays to zero as L —>• oo, and hence 
E[e Ujv ^] ~ (l — F (-ny 1 ) ) r ( x + «) as N -> oo. 

By a similar argument one obtains that 

lEfuTv^e^^] ~ (l — F (—My 1 )) (T(l + a) — T(a + 2)) as N —>• oo. 
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The formula unA'(un) — A (un) = In N yields 

(i - F (-»«)) N 


i 


as N 


e a aT(a) 

We now use (11.7p and Karamata’s representation to conclude that 


oo. 


u 


-1 

N 


lim 

N —>-oo Cl]\r 


1 


oT(a) 


1/a 


and hence as N —» oo 


A~-~ —— 

mat e 


ctT(a) 


l/a 


Oat, 


which proves the statement. The second claim is a direct consequence of (13.21) . 

If / 0, we can simply translate the by x$, so the hypothesis of Proposition 
In the next corollary, we prove the upper bound in Theorem 11,11 


Corollary 3.3. Assume that the hypothesis of Theorem [777] hold. Then, as N 


oo 


□ 

hold. 


limsup('c,v(^) — x^a^f < — c a . 

N—foo 

Proof. In the case x £ = 0, it is a straightforward consequence of Proposition 13.21 and (13.21) . If 
X£ 7 ^ 0 , it suffices to translate the variables fij by x^. □ 


4 Lower bound 


In this Section, we show that for every e > 0 there exists a N$ such that V/V > /Vo 


(Ujy(0 -xg) 

qn 


> -C a - £, 


(4.1) 


which proves the lower bound in Theorem 11.11 Throughout this section, we fix an arbitrary 
e > 0, and we assume that = 0 without loss of generality. To prove (14.11) . we construct a 
process x(t) € that bounds X{t) from below, hence 


max (x(tf) < max (X(f)). 

Then, in Subsection EH we check that the process x(t) is a M-BRW, and we show that for M 
large enough and the appropriate offspring distribution (see Subsection 14.2(1 , 

lim inf - max (x(t)) > —(c a + £)ajy a.s. 
t —>oo t y ' 

which implies (14.11) for x^ = 0. The general case is obtained by a simple affine transformation. 
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4.1 Coupling with a Branching Random Walk 

We construct x(t) inductively as follows: let M, K € N, the appropriate values for K and M will 
be chosen later on, and assume that N > KM (in fact, we will take KM negligible compared 
to N). For t. = 0, we define 

Xi{ 0) = X ai (0), 

with <7; = <7j(0). Assuming that the process x(-) has been constructed up to time t £ N, the 
vector x(t. + 1) € is obtained according to the inductive rule below. 

1. Branching step: Every particle Xi(t ) is replaced by K new particles (reproductive law), 
whose positions are defined by a point process (xj(f)) translated by Xi(t). 

The point processes (xi{t )); 1 <i<M) are also constructed according to an induc¬ 

tive rule, that we describe: 

• For i = 1, let 71 := {1,..., N — KM} and denote by 

4iw(* +1) > Clb* + !)>•••> +!). 

the K largest values among {^ai(t),j (t + 1); 3 £ 7i}. Let, also, 

Ji = Ji(t + 1) 

be the set of their indices, that is, 

we will keep track of these labels. Note that the indices jj '■* = 1 <l< K 

are random. Then, i(t)) is the point process 

(t)):= J2 1)}’ 

JC (1+1) 

obtained by the K largest values in 7i and the descendants of xi(t) are at the 
positions: 

xi (t) + (t + 1) for 1 < l < K. 

• Assume that we have constructed l<j<i — l) and that we have kept 

track of the sets J\,... appearing in the respective constructions. The sets 

Jj = Jj(t + 1) C {1,... ,N} are random and disjoint. Then, given J\ U ... U Ji-i, 
we choose 

Ti = Ti(t + l) c {l,..., N} \ (Ji u ... u Ji_ i) 

according to a deterministic rule. For example, one can choose the N — MK first 
elements (in the usual order of N) in {1,... , N} \ ( j\ U ... U 3%-\) ■ By construction, 
Ti is a random set of {1,..., N} satisfying the property 

T i nj 1 = Q = T i nj 2 = ... = T i n Ji-i. 
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Let 


t<\-Ti) 


(t+i) >i 


(2=75) 

<H (t) 




be the K largest values among {Co*(t),j(* + !); 3 € %}, and Ji(t + 1) = {jj' 1 ,..., j$} 
be the set of their indices, that is, 


t(i:75) 

i(t ) 


(* + 1 ) “ + 


Then, C^{xi (£)) is the point process formed by the these K points. 


We end up the branching step with KM new particles. 


2. Selection: We select the M rightmost particles among the KM obtained in the branching 
step. 


3. Ordering: We reorder the M selected particles to obtain the vector x(t + 1). 


In the next two lemmas, 
that the point process 


we show that x(t) -< X(t) and that are i.i.d. which implies 

M 

i =1 


has the distribution of the point process obtained from a M-BRW 

First, we prove that x(t) bounds X(t) from below. We bring to the reader’s attention that the 
next lemma is a direct corollary of the construction of x(t) and it holds without any assumption 
on the family {^ (s); 1 < i, j < N, s > 1}. 


Lemma 4.1. For N > MK, let x(t) be the branching/selection process constructed as above. 
Then, x(t) bounds X(t) from below. 


Proof. It is immediate that x(0) -< X(0), hence assume that x(t) -< X(t). Before the selection 
step, there are MK points at the positions 


Xi(t ) + £ . . .(i) (t + 1) 1 < l < K and 1 < i < M. 

By the construction of x(-) the jf' 1 are all distinct. Since Xi(t) < X a .^(t), we have that 

Xi (*) + + ^ - X j^ + 

Hence, the point process obtained from the branching step bounds X(t + 1) from below, as a 
consequence, after the selection step x(t + 1) ~< X(t + 1), proving the statement. □ 

Now, we prove that the point processes are i.i.d. Lemma 14.21 holds under the unique 

assumption that the family {^ (s); 1 < i, j < N, s > 1} is i.i.d. 
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Lemma 4.2. Assume that N > KM, and let { %(t);t 6 N;i = be the set of indices 

obtained in the above construction. For t > 0 denote by ’E[Xi(t)) the point process 


[Xi 


wb= E 


j’S7i(t+l) 




then, |S(xj(t)); 1 <i< M; teN} are i.i.d. 

In particular, the family of point processes (xi (t)); t € N and i = 1,... , M} is also 

i.i.d. 


Proof. Note that the families of random variables 

{<r(s); 0 < s < t}, (7i(s);0 < s < t; 1 < i < M}, {H(xj(s)); 1 < s < t — 1; 1 < i < M}, 

are ^-measurable with Tt from (12.9|i . By assumption, a { fij (f +1); 1 < i, j < TV} is independent 
from J-j, then, by successive conditioning, one easily checks that conditionally on Tt the vector 

(Co-+ 1); i = 1, • • • , M and j € %) 

is distributed according to a Af x (iV — KM) vector, whose entries are i.i.d. copies of , which 
implies the independence from 75- Moreover, the conditional independence of the + 1) 

yields that (E(xj(t)); i = 1,..., A/) are also independent, which proves the first claim. 

The second claim is an immediate consequence of the first part of the lemma. □ 

Finally, we focus on the asymptotic distribution of C^ K \-) after suitable rescaling, 

■■= E Wr (4 - 2 > 

zeC( K ) (xi(t)) 

for a at given by (11.61) . With some abuse of notation, we denote by 0 >( N ’ K ) the common distri¬ 
bution of these point processes. 

Proposition 4.3. Assume that the hypothesis of Theorem [777] hold with x^ = 0 and that M and 
K are fixed. Then, as N —>• oo, 

£*( N ' K ) —■> 0 > (K) in law, 

with defined in Corollary \2.2 1) with (5 = a — 1 and C = a. Moreover, for every t > 0 the 

moment convergence 


lim E 

min 

£' 

= E 

min ,5^*) 


TV—>-00 







also holds. 

Proof. It suffices to prove the convergence for IP^ N,K ){x\ (t)). Since is in the domain of 
attraction of $ a and x £ = 0, for every z > 0, 

IP (C<n(t)j(*) > -™ n ) ~ ^"1P > ~«iv) ~ as IV -)• oo. (4.3) 
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It is a classical result of extreme value theory [22] that, as N —>• oo, the point process 


gt>(N,K) 


N-KM 


law 


3 = 1 


converges in distribution to a PPP with intensity measure \z\ a ^ 1 a\i z< ^dz, as claimed. A 
necessary and sufficient condition for the convergence of the £th moment is that the r.v. 
has itself finite tth moment, which is a consequence of |(H1)[ Proposition 2.1 in [22] proves this 
statement for the maxima of i.i.d. random variables in the domain of attraction of 'P Q . 

Now, a line-by-line adaptation of Proposition 2.1 in [22] yields the last claim. □ 

A straightforward consequence of Proposition 14.31 and the two previous lemmas is that 


K. 


(N) 


(t) := 


M 

2—1 




converges in distribution to the point process obtained from a M- BRWY^^)) at time t, more¬ 
over, 

max (^ (7V) (*)) < \ ^ ■ 

We will prove in Subsection 14.21 that if one chooses K and M large enough (depending only on 
e and the distribution ), then for N larger than some No > 0, 

liminf ^ max (^ N Ht > — c a — e a.s. (4.4) 

which proves the lower bound (14.11) . 


4.2 Uniform lower bound for the velocities 

In this subsection, we prove the lower bound (I4.4p . which concludes the proof of Theorem 11.11 
The proof is divided in two main steps. 

In the first one, we focus on the BRWs defined by &( N ’ K ) _ yy e p r0 ve that if N is suffi¬ 
ciently large, with positive probability there exists more than M vertices w in generation n (see 
Subsection 14.2.11 for its definition), such that 

position (wt) > — (c a + e/2)t \/wt € [e,u;]. 

In the second step, we use this result to obtain the uniform lower bound (14.41) for the Af-BRWs. 

Most of the arguments presented here have already been used by Berard and Gouere [3]. In 
our case, though, we deal with a family of BRWs indexed by N. whereas in [3] they compute 
the velocity for a unique M-BRW. 
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4.2.1 First step 


Let ^^ N ’ K ) be the distribution defined by (14.21) and denote the distribution of a point 

process obtained from the K largest points of a PPP with intensity measure \z\ a ~ l a\^ z< ^dz. 
Then, BRWf^^’^) and BRWf^^l generate the same Galton-Watson tree, in which every 
individual has a constant number K of offspring, denote by Tk this tree. We will construct 
these BRWs on a same probability space. 

Let {^ N,K \w)]w £ Ta-} be i.i.d. copies of {&( K \w)] w £ T/<} be i.i.d. copies 

of and (0,7-*, P) be a probability space where those families of r.v. are defined. Since 

cp(N,i\) con verges to in distribution (see Proposition 14.31) . we can and we will assume that 

the stronger a.s. convergences 

lim &^ N ' K \w) = &^ K \w) a.s. (4-5) 

N^-oo 


hold for all w £ Tk : which implies the point-to-point convergence 


lim 

TV—>• oc 


@>( N ' K \w) - &>W(w) 


= 0 


a.s. 


where || • || is the distance defined in (EH). Note that we have not lost in generality, since we 
can always construct a probability space (fl, T, P), for which the a.s. convergence holds. 

Under these hypothesis, the construction goes as follows. Each individual w £ T k has K 
offspring, that we label according to some deterministic order. Let be its ith children, then, 
its position rj( N \ww^) and r](°°\ww^) in the BRW(l 3 2i ^ Ar ’ A ^ an d BRW(7^ a ^) are given by 

(ww^) = r]( N \w) + ZP\ N ' K \w) and rf°°^ (ww^ 1 ^ = rf°°\w) + K \w ), 

where ^ N,K \w) and &[ K \w) denote the ith largest point in £?^ N ' K \w) and ^ K \w) re¬ 
spectively. This construction couples the BRWs and for w £ T k fixed 

lim rf N \w) = rf°°\w) a.s. 

7V->oo 


A direct calculation shows that satisfies (I2.4H2.6D . which implies the existence of the 

asymptotic velocity 7 (t 5 2> ( A '), with 7 given by (12.71) . Lemma \2 . 1 1 with C = a and f3 = a — 1 > — 1 
yields 

lim 7 ( 7 ^^) = —— ^ 


K —>00 


r(a)a 

Let 6 = e/12, then there exists Kq such that V/L > I\q 

7^^(R)) > _ Cq _ § 


= ~C n 


(4.6) 


Fix K for which (14.61) holds; we bring to the reader’s attention that, as (HI). (H2) hold, £P( N ’ K ) 
also satisfies {2.1$ , \2. 51) . Moreover, a simple calculation shows that 7 (^ 2 'L V ’ K )) converges to 
7 (^W) 

as N —>• 00 . 

We now prove that with positive probability there exists more than M individuals w £ T in 
generation n such that 

(u’t) > —c a t — 65, for every -wt £ [e, fi)]. 
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As it will become clearer in the sequel, we take n of the form n = sw + m, with 


'In AT 

+ 1 and m = 

(|-R| - c a - 65) s M 

In p 


3 5 


(4.7) 


The constants ip > 1 and R < —c a — 65 < 0 are given by Lemma [4.41 and formula (14.81) below and 
depend only on the distribution ^ K \ Although M may be very large, it will be kept constant 
throughout this section (while N —> oo), hence sm and m are also constants. 

First, we obtain a lower bound for the probability of the set 


1 3w E T k in generation m such that r/^ N ^{wt) > (— c a — 3 5)t; Vwt E [e,n;]|. 


Denote by A m $ the set 


A m ,5 ■ 


| £p>( N J<) ( u /) _ < <5;Vu/ E T k such that | w'\ < m| 


then, for m E N and 5 fixed, one obtains from ( 14 . 51 ) that P(A mj , 5 ) —>• 1 as N —>• oo. Since 
7 (<^ (A) ) > c a — 6 we have the following set inclusions 


{3u> E T k such that |tc| = m and 77 ^ (wt) > (— c a — 3<5)t; \/wt E [e,rc]} 

D E T k such that |w| = m and (wt) > (— c a — 3 <5)t; Vuy E [e,u;]} n A m j 
D {die E T k such that |u;| = m and > (—c Q — 2S)t; 'iwt E [e,to]} H A m 

D E T/v such that |tc| = m and r)(°°^(wt) > — 5)t; \/u’t E [e,re]} n A m ^ 

D (BRW(^ (a) ) has an infinite ray lying above the line of slope — 5)} n A m ^, 


which yields the lower bound 

lim inf P (die E T k such that liol = m and rj^ (wt) > (— c a — 3<5)t; Vwt E [e, wi) > p(oo,5). 
N-> oo V v ' / 


From Theorem 12.31 p(oo,5) > 0 is a constant depending only on ^ K \ Then, there exists 
Nm E N depending only on rri (and hence, on M) such that \/N > Nm 

P^die E T k : |to| = m,rf^ (wt) > (—c Q — 3 S ) t , Mw t E [e,tc]^ > ^ . 

Now, we choose R and <p in (14.71) . Since (] — oo, 0[) = K, one can take R < —c a — 65 < 0 
such that 

f(&>( k) [R,0) > 2) > (4.8) 

Using the convergence in distribution, there exists a N' > 0 such that for N > N' 

P (&( N ' K \R, 0) > 2 ) > jj. 

Without loss of generality, we can and we will assume that Nm > N'. The Galton-Watson tree 
whose offspring distribution is 

p\ N) = P ^ N ’ K ) [R,0) = i), i = 0,1,..., (4.9) 

has mean offspring larger than 4/3. It is supercritical, and the following well-known result holds. 
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Lemma 4.4 ([2] Theorem 2 Section 6 Chapter 1). Let Mt denote the population size of a 
supercritical Galton- Watson process with square integrable offspring distribution (started with 
one individual). Then, there exists r > 0 and <p > 1 such that for all t > 0 

P (M t > </) > r. 

Let Mj N> denote the population size of the Galton-Watson processes defined by (p\ N ^)i= 
Using a simple coupling argument and Lemma [4.41 we can find a ip > 1 and r > 0 not depending 
oniV> Nm such that for all t > 1 

P {M ( t N) > tp*) > r. 

With m and sm from (|4.7I) . we have that 

> M ; with probability at least r > 0, 

and that 


(— c a — 3 5)m + R(t — m) > (— c a — 6 5)t; for every m < t <m + sm- 
Let ww' G Tk be a vertex in generation n, with the following properties: |u;| = m , 

7] (N) (wt) > (~c a - 3 5)t, Mw t G [e,w], 
w' is in generation sm in the Tj< sub-tree descending from w and 

v {N) («4+i) - v {N) K) > R . Vw's e IK w'l 

Then, by a simple calculation one can conclude that the path [e, ww'} C Tk has always lain 
above the line of slope — c a — 65. For N > Nm, a conditioning argument yields the lower bound 
for the probabilities 

P |tc G T k', \w\ = n and rf' N \ wt ) > — [c. a + t Vwt G [e,u)]| > M^j 
> P^3rc G T k such that rf N \wt) > (— c a — 3 5)t; Vwt G [e, w} and M 

>r P (°°’ S ) 

2 ’ 

in the second equation, w is a vertex in generation m and is the population size of the 

Galton-Watson process generated by the descendants of w for which 

v {N) K+i) - p {N) («4) > R ^ w 's e K w'l- 

In particular, we have just proved the following proposition. 

Proposition 4.5. Let w€Tk) the BRW defined by the point processes LpNN ). Given 

£ > 0 let R be given by 0, r and <p as in Lemma \4-41 Then, take sm and m as in (B- 
Then, with n = m + sm, we can find some Nm (depending only on M) such that if N > Nm 

P jw G T K ] |w| = n and r/ (N \w t ) > - (c a + t; \/w t G [e,w]| > M^j > 
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4.2.2 Second step: uniform lower bound for the speed 


In this step, we obtain a uniform lower bound for the speed of a^xit), which is simply a 
M-BRW(^ 7V,A )). Let ^ N \t) be the point process associated to a^x(t) 

M 

2—1 

and 7 m(^ (W ) be the asymptotic velocity of the M-BRW(^ iV ’' K '}), 

^lim t~ 1 min (^- N \t)^ = ^lim f _1 max (^ N \t)^ = 7 'm{^ N,K ^) a.s. 

Then, we will prove that for e > 0 and K given by (14.61) the inequality 


liminf (liminfy ’m(^ N ’ K ^)) > — (c a + e) 

iW-Kx> \ 7 V-kx> ) 

holds, which finishes the proof of Theorem 11.11 Following the strategy of [3], we construct a 
third point process W(t) = W^ N \t) that bound ^ N \t) from below. This new point process 
evolves like x.^ N \t) up to a certain random time r,, i G N, from which we shift the position of 
all particles to the minimal position, and start W ( t ) afresh. 

Let n = m + sm , where sm and m are given by (14.71) . We will construct the process W(t) 
and the stopping times 0 = tq < n < ... together 

7~i := inf jl < s < n; min > (— c a — e/2)s| , 

where inf{0} = n. Then, ri < n is a stopping time with respect to the filtration J-f. For 
0 < t < t\ let 

and mi := min (W(n)), then at the time step ri —>• t\ + 1 we shift all particles W t to mi and 
continue the construction up to T 2 according to the induction step. 

Inductive step: assume that T\ <... <ti and W(t) for f <77 are defined. Then, for r/+l<t< 
T] + \ (we will define 77+1 below), W(t) is the point process of a M- BRW(^ iV ’ A )) starting from 

mi := min (W (t/)) . 

At each time step t —>• t + 1 the individuals (Wj(t )) i=1 M give birth to K new individuals, 
whose positions are determined by independent point process ^^ N,K \xi{t))\i = 1 ,...,M), 
and die immediately afterwards. We assume that the point process defining ^ N -* and W(-) are 
the same. Moreover, we will also assume that the indices are organized in order to couple 7 ^ 
by W(-). We then select the M rightmost particles to form the next generation. 

The process evolves as above up to 

n+ 1 := inf {ti T 1 < s < n + n; min (W (s)) - m t > (~c a - e/2)s], 

where we shift the positions of the M particles to mi+ 1 , the minimum of the positions. It is 
immediate from the construction of IT (■) that 

W(t)^z {N \t) VteN. 
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For l > 1, the processes (W(t) — mf, t € [r; + l,r; +1 ]) and the random variables r / +1 — 77 are 
i.i.d. In the sequel, we use the notation r := n, then by the law of large numbers 

^lim y min ^^( 77 )^ = 7 Af (■^ JV ’' K ^) 1 E[t] a.s. 

From the construction of W(-) and the renewal theorem we also obtain that 

lim inf — min (ti) \ > lim inf — min (W ( 77 )) = E [min (IF(r))] a.s. 


which implies that 

■mi*™) > E[miD E { T v [ - ' (T>)] . 

With B = { min (lF(r)) < (— c a — e/2)r}, we write 

min (IF(r)) > (-c Q — e/2)rl B c + min {W(n))l B 

= (—c Q — e/2)r + (c Q + e/2)t1 b + min (!F(n))lB. 


(4.10) 


Taking expected value we get 


E[min(lF(r))] > (-c a - e/2)E[r] + E[min (W{n))l B \ ■ (4.11) 


Let min y&y N ' K > (lF,(t)) J be the smallest point of the point process generated by IFj(t) before 
the selection step, it has the law of the K th maxima of a IV — KM sample of Since < 0, 
one gets the lower bound 


n M 

min (IF(n)) > min (iFj(I))^ , 

£=0 z=l 


which implies that 


E[min (lF(n))ls] >— (n + 1)ME 


min 



The Cauchy-Schwarz inequality yields 


E[minIF(?r)lB] > — (n + 1)ME 


min 


0>{N,K) 


21 V 2 

E(R) 1/2 . 


By Proposition 14.31 the second moment of min converges as N —>• 00 to a finite 

constant. Hence, there exists a constant c, depending only on ^, such that 

E[minIF(n)lB] > —c(n + l)MP(R) 1 / 2 . 

Finally, the probability of I? can be estimated using Proposition 14.51 The evolution of different 
individuals in the Ai-BRW is not independent. Yet, a M-BRW can be coupled with M inde¬ 
pendent BRWs, see Section 3.3 in [3], so that the event “the minimum of the M-BRW always 
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lies below the line of slope — c a — e/2” implies that none of the M independent BRWs has more 
than M vertices in generation n that have always stayed above this line, hence 


P(B) <11- 


M 


From the definition of n, if M is large enough c(n + 1 )M < M 2 and 


lim sup M ( 1 — 

M—>• oo 


p{ oo,e/12) \ M/2 


= 0 . 


Then, choosing M properly (note that it depends only on e and £ t j but not on N), one gets 

E [min W (n) 1b] > — 

Combined with (14.10114.1111 , this ends the proof. 

Acknowledgements: The authors thank Bernard Derrida for stimulating conversations on 
this model. 
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